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We consider some of the methods that can be used to reveal the general features of how wave 
functions evolve with time in the harmonic oscillator. We first review the periodicity properties over 
each multiple of a quarter of the classical period of oscillation. Then we show that any wave function 
can be simply transformed so that its centroid, dehned by the expectation values of position and 
momentum, remains at rest at the center of the oscillator. This implies that we need only consider 
the evolution of this restricted class of wave functions; the evolution of all others can be reduced to 
these. The evolution of the spread in position Aj, and momentum Ap throws light on energy and 
uncertainty and on squeezed and coherent states. Finally we show that any wave function can be 
transformed so that Aj, and Ap do not change with time and that the evolution of all wave functions 
can easily be found from the evolution of those at rest at the origin with unchanging Aj, and Ap. 


I. INTRODUCTION 

The simple harmonic oscillator is studied in every 
course in elementary quantum mechanics because of its 
relative mathematical simplicity and its usefulness as a 
model of natural processes: oscillators are ubiquitous in 
nature and harmonic traps are used in many labs. Yet 
little attention is given to the time-evolution of oscillator 
wave functions, beyond the fact that the evolution is pe¬ 
riodic. Many texts do derive the propagator, but using 
it requires an integration that can be carried out exactly 
only for a small class of wave functions.^ 

Here we will consider some of the simple methods that 
can be used to understand the general features of how os¬ 
cillator wave functions evolve. We start with the simple 
relations between the evolved and the initial wave func¬ 
tions after each multiple of ^T, where T is the period of 
the equivalent classical oscillator. Then the evolution of 
some moments of the wave functions will be determined. 
The first moments, that is the expectation values of posi¬ 
tion and momentum, give Ehrenfest’s result that (x) and 
(p) follow a classical trajectory. A powerful result is that 
the evolution of oscillator wave functions is essentially 
independent of the first moments. Consequently we need 
consider only states that remain centered on the origin 
(taken to be the point where the oscillator restoring force 
is zero). All other states can be reduced to these by a 
simple transformation. 

The second moments give a rich set of results about un¬ 
certainties and energy. Squeezed states emerge as Gaus¬ 
sian states at the extreme edge of a strong form of the 
uncertainty relations, an inequality involving the second 
moments. All Gaussian states can be produced by chang¬ 
ing the spatial scale of (that is squeezing) the ground 
state, and applying the transformations discussed in the 
previous paragraph. Similarly, any state can be pro¬ 
duced, apart from a phase, by changing the scale of its 
‘stable’ form whose second moments remain constant as 
it evolves. The evolution of any state can be obtained 
from that of its stable form. 

A note on the use of dimensionless variables: the oscil¬ 
lator has its own natural length scale a = and 


the theory becomes much simpler in terms of the dimen¬ 
sionless variables xja, apjfi and uit. The oscillator has a 
symmetry between position and momentum (essentially 
because each appears only as a quadratic term in the 
Hamiltonian) and this is somewhat obscured when us¬ 
ing the usual physical quantities x and p. Almost every 
equation in this paper would be simpler, and the gen¬ 
eral arguments clearer, if expressed in the dimensionless 
variables. The downside would be that it is harder for 
students to relate the oscillator theory to that of other 
systems if the discussion is in terms of special oscillator 
variables. For that reason physical variables will be used, 
except for Eq. ^ of Section HI and all of Appendix B, 
which involve Fourier transforms. 


II. PERIODICITY 

For the oscillator with mass m and angular frequency 
w, the Hamiltonian is 

H = —-I- (1) 

2m 2 

and the energy levels, that is the eigenvalues of H, are 
En = hu}{n -\- ^) with n = 0,1, 2,.... If a wave function 
'>f{x) is expressed as a sum over the energy eigenfunctions 
ipnix) at time t = 0, so that if^x) = Cn'4’nix), then 

the evolved wave function at time t is 

ip{x, t) = ^ exp[-iLut{n + -)]c„tpn{x), (2) 

n=0 

where i = a/—! and the c„ are generally complex. 

Evolution over one period. If t = 27r/a;, the factor 
exp[—iojt{n+^)] becomes —1, so that if{x,T) = —ijj{x,Q) 
where T = 27r/a;. Thus, every wave function returns 
exactly to its original form after each period T, except 
that the wave function changes sign. This period T is also 
the period of the corresponding classical oscillator. This 
periodicity is a direct consequence of the equal spacing 
of the energy eigenvalues. 
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Evolution over half a period. At time t = tt/uj, 
exp[—iojt{n + |)] becomes — 1 (—1)". Since the energy 
eigenfunctions are alternately even and odd, ' 0 „(—x) = 
(—l)"^„(x), it follows from Eq. § that 

tjj{x,T/2) = (3) 

This holds for any initial time and therefore 

tpiXjt + T/2) =—i-tjj{—x,t). (4) 

Evolution through the second quarter period. For 

an initially real wave function, the evolution for all time 
can be deduced from that over the first quarter pe¬ 
riod. Inserting time t' = T/2 — t into Eq. (1^ leads 
to '4>{x,t') = —l)”c„ exp[ia;<(n -|- Since, 

4’n{—x) = (—!)"''!/'»(a;), and the c„ are real, we have 

ip{x,T/2 — t) =—lip*{—x,t) if0) real. (5) 

Thus, after one quarter period, \ip(x,t)\ retraces the val¬ 
ues it took during the first quarter, but it is spatially 
reversed and the real part takes the negative of the val¬ 
ues that the imaginary part took. Only the evolution 
over one quarter period need be calculated and the rest 
can be deduced. An example is shown in Fig. 1 with 
more detail in Section IX. 


T t=T/4 



FIG. 1. 

The evolution over a half period of a superposition tp of two 
Gaussian wave packets, each with the same width as the 
ground state. The larger packet starts at rest with its center 
at X = 20a while the smaller, with amplitude 0.4 of that of 
the larger, starts at 17a and also at rest. The solid curves 
show |i/)| for five times over a half period while the dashed 
curves show the magnitude of each Gaussian if the other 
were not present. The thick horizontal lines show the spread 
in position of the superposition. The initial wave 
function is real and therefore the graphs for negative x are 
the mirror image those for positive x, as in Eq. ([^. 

Evolution over a quarter period. At t = 7 r/ 2 a;, 
exp[—itcfn] becomes (—^)" and this can be compensated 
by taking the Fourier transform of the eigenfunctions, as 
shown in Appendix B. Consequently every wave function 
changes essentially to its corresponding momentum wave 
function over f T. This result emerges directly from the 
propagator and is discussed in the next section. 


III. THE PROPAGATOR 

In principle, the problem of the evolution of wave func¬ 
tions is solved by finding the propagator K{x,x',t) such 
that, for any initial wave function ip{x), its evolution at 
time t is 

/ OO 

K{x,x',t)'tfj{x')dx'. ( 6 ) 

-OO 

Derivations of the propagator for the oscillator have been 
given in many texts [for example. Ref. [H Section 15.2] 
and paper^ and a derivation is given in Appendix A. 
The result is 

K{x,x',t)= . , L^^ X (7) 

OLy2Txi sin cut 

exp] ^ „ -((x^ -I- x'^) cos wt — 2xx')l, 

2a‘‘sin Lvt 

where a is our length scale The propagator 

contains t only as sinwt or coswt, and is therefore peri¬ 
odic with period T = 27r/a;, but there is an ambiguity 
i n phase d due not only to the ambiguity in the sign of 
Vsin uit but also to the singularity in the exponent when¬ 
ever sinwt = 0. This ambiguity is resolved in Appendix 
A; the result is that an extra factor of —t must be in¬ 
cluded each time sin cot goes through zero. 

A simple result that can readily be derived from the 
propagator is that every Gaussian wave function remains 
Gaussian as it evolves. 

The propagator also gives a simple result after any 
quarter period of evolution: at time t = T /4, 

1 tt' 

'ijj{x,T/A) = — 7 ;=/ exp[-i—]'tp{x')dx'. ( 8 ) 

aV27r* J_oo « 

Because of the oscillator’s symmetry between position 
and momentum, it helps to use the dimensionless vari¬ 
ables ^ = x/a and p = ap/h when taking Fourier trans¬ 
forms of wave functions. Write'0(x) = 4 '(x/q;) and define 
its Fourier transform as 

^(P) = f exp(-ip^)4'(^)(i^. (9) 

V J — OO 

Then ip{x,T/A) = exp(—^ 7 r/ 4 ) <i)(x/a). Thus, any state 
evolves over a quarter period into a state whose spatial 
wave function equals the Fourier transform of the initial 
wave function (apart from a phase change of —it/A). 

Thus the propagator can be used to draw out some gen¬ 
eral features of the evolution. For application to specific 
initial wave functions, however, the integral in Eq. is 
difficult to calculate numerically because the integrand 
has rapid oscillations in some regions. Other methods 
are used for numerical work (such as using an expan¬ 
sion of the wave function over the energy eigenfunctions). 
Furthermore, the integral can be evaluated analytically 
only for very few forms of initial wave function. Other 
than for an initial Gaussian (or higher energy eigenfunc¬ 
tions) those that can be evaluated give results in terms 
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of special functions from which it is difficult to deduce 
the general features. For example, an initial square wave 
function can be evaluated in terms of a combination of 
error-functions with complex arguments.® Note that the 
square wave function is too singular to be subject to most 
of the analysis below because its energy is infinite. 


IV. DYNAMICS OF THE OSCILLATOR 


We will use the notatiorP 
dA 


DfA = 


dt h 


[H,A]. 


( 10 ) 


Then DtA is the ‘total time-derivative’ of the operator A, 
which adds to the partial derivative a term that takes ac¬ 
count of the evolution of the state under the Hamiltonian 
H in such a way that 


dt{A) = {DtA) 


( 11 ) 


for any state and any operator A, where dt stands for the 
usual time-derivative d/dt. It is easy to see that 


Dt{AB) = {DtA)B + A(DtB), 


( 12 ) 


similar to differentiation of a function. 

For the oscillator Hamiltonian in Eq. Q, it is simple 
to calculate that 


Dtx = p/m, Dtp = —moj^x, 


(13) 


similar to the equations of motion for the classical os¬ 
cillator. These equations will be used to calculate how 
various expectation values evolve for the oscillator. 

Motion of the centroid. We will refer to the expecta¬ 
tion value of position as the centroid of the wave function. 
From Eqs. © and ( [T^ 


dt{x) = {p)/m, dt{p) = -mJ^{x). 


(14) 


These are the usual classical equations for the oscilla¬ 
tor with the solution {x) = A cos uj{t — tp) and (p) = 
—TOW A sin a; (t — to), where to is the time of maximum 
displacement A. This is Ehrenfest’s result that the cen¬ 
troid follows classical sinusoidal motion. 


ip{x), we can find the initial position Xq = (x)^ and mo¬ 
mentum po = {p)tp of its centroid. Then the complete tra¬ 
jectory of the centroid follows from the classical equations 
of motion: dtX = p/m and dtp = —muAx. The centroid 
position is x{t) = xq cos uit + (po / tow) sin wt and momen¬ 
tum is p{t) = Po cos wt — mujxo sinwt. We can produce 
from ip{x) a wave function (/{x) = exp{—ipQx/h)'ip{x+Xo) 
whose centroid sits at the origin with no momentum. If 
(/{x, t) is the evolution of this initial wave function, and 

i/ix, t) = exp[^p(t)(a: - |)]((i(a; - x{t),t), (15) 


then, with ^ = a: — a;(t), it is not difficult to show that 

(16) 


h? d'^ I 2 2 ^ 1 / / \ 

-I- -muj X — ih—\ih(x,t) 
^ 2m dx^ 2 ’ 


=exp[ifi(()(,T - f 

That is, if 4>{x, t) satisfies Schrodinger’s equation, then 
so does '4>{x, t). 

Displaced eigenstates. As a simple example, consider 
the displaced ground state tpix) = exp[—(a: — a)'^/2a^], 
where a is any real constant. The centroid will follow 
x{t) = acoswt, p{t) = —TOwasinwt. With ^ = a: + o 
we have (^(^,0) = exp[—^^/2a^], which is the ground 
state whose evolution is (l){^,t) = exp[—^^/2a^ — iojt/2]. 
Therefore the evolution of ijj is, from Eq. (15), 


ip{x, t) = exp[i0 — {x — acos wt)^/2Q;^], (17) 

where 6 = —asinwt (a: — |acoswt)/a^ — ^wt. 

The same is true of the excited states. As shown in ev¬ 
ery textbook, the eigenstate has the form tpn{x,t) = 
exp{—x‘^/2a‘^)Hn{x/a)exp[—i{n+ ^)wt] where Hn{z) is 
the Hermite polynomial. It then follows that there 
is this much larger class of displaced excited states, each 
defined by n and the initial position ao and momentum 
Po of the centroid, that have the evolution 

4'„(a,t) = exp[i6»„ - ^^^/a^]i7„(f/a), (18) 

where ^ = x — x{t) and 9n = [p{t)^ + ^p{t)x{t) — Ent]/h. 
Here = huj{n+A) and x{t), p{t) give the classical evo¬ 
lution of the centroid, x{t) = xocosuit -I- (po/TOw)sinwf 
and p{t) = pocoswf — TOwaosinwt. For each of these 
states the magnitude |'I'„(a — a, t)| does not change with 
time, while the centroid a follows a classical oscillation. 


V. EVOLUTION IS INDEPENDENT OF THE 
CENTROID 


VI. EVOLUTION OF THE SECOND MOMENTS 


Here we will show that the evolution of any wave func¬ 
tion is independent of the motion of the centroid, except 
for a simple effect on the phase of the wave function. This 
greatly simplifies the problem of finding the evolution of 
an arbitrary wave function: one need only consider the 
evolution of wave functions whose centroid sits at the ori¬ 
gin and does not move. Given any initial wave function 


It is worthwhile to study the second moments of the 
wave function, 


= (a^) - (a)^ 

= if) - (P)^ 

Axp = X (pi + ip) - (i) (p), 


(19) 

( 20 ) 


( 21 ) 




4 


because A^; and Ap are measures of the spread in position 
and momentum of the state, while A^p is a measure of the 
correlation between position and momentum, but also 
because these moments enable insight into the energy 
and uncertainty properties of the state. We can use the 
operator equations (12) and © , to derive 


Dtx^ = {px + xp)lm, (22) 

Dtp^ = —muj^{px + xp), (23) 

Dtipx + xp) = 2{p^ jm — muj'^x^), (24) 


and there is a similar set of equations for (x^), {p^) and 
{px + xp), as well as for (x)^, (p)^ and 2(p)(x). Hence 

dtAl = 2A,p/m, (25) 

dtAp = -2mu}‘^A^p, (26) 

dtAj-p = Ap/m - muP'Al. (27) 

It follows that d^Axp = —dw^A^jp and therefore A^p 

varies sinusoidally with twice the frequency of the os¬ 
cillator. We can take 


Axp = hAsm2u!{t — to), (28) 


where tg is a time when A^jp is zero. Then, integrating 
Eqns.(25) and (26), 

A^ = a^[e — Hcos2a;(t — fo)] (29) 

Ap = (h/a)^[e + A cos 2w(t — to)]• 


The constant e must be the same in each of these equa¬ 
tions to satisfy Eq. ^7^ . These basic properties are illus¬ 
trated in Fig. 2 showing tPpjo?, a^Ap/h^ and Aa;p/h over 
two periods of the oscillator, for the sum of two Gaus- 
sians shown in Fig. 1. More detail is given in Section 

m 

The question of how to calculate the constants e, A 
and to in terms of A^, Ap and A^jp will be taken up at 
the end of the next section. 

In summary, for any wave function, the evolution of the 
second moments can be simply calculated from their val¬ 
ues at any initial time. The spreads A^ and Ap oscillate 
sinusoidally about fixed values at twice the frequency of 
the oscillator, always differing in phase by 180 deg. The 
correlation A^jp oscillates about zero with phase differing 
by 90 deg from the others. Although A^, and Ap can be 
interpreted as spreads about the centroids, their evolu¬ 
tion does not depend on the motion of the centroid. 


VII. ENERGY AND UNCERTAINTY 

We define the quantal energy to be Eg = (H) — Ec 
where Ec = (pf' 12'mA muP'(x)^ 12 is the classical energy 
of the centroid. Thus 

^ (30) 

The constant e = Eg/tux) is then the quantal energy in 
units of hui. Both Eg and Ec remain constant. 




FIG. 2. 

The oscillations of the first and second moments of the 
system described in Fig. 1. The top graph shows the 
position (x) and momentum (p) of the centroid, each scaled 
to be dimensionless, plotted against the time in units of the 
period T = 2n/uj of the oscillator. The lower graph shows 
A^, Ap and Aa;p, each scaled to be dimensionless. The 
amplitnde of the oscillations in each of these three cnrves is 
A. The vertical scale for the centroid (top graph) bears no 
relation to that for the second moments, although the time 
is aligned. In general, the phase of the oscillation of the 
centroid need not align with that of the second moments. 


Although the initial values of all three second-order 
moments are needed to calculate the evolution of any of 
them, the combination 


K^ = {AlAl-Al)in^ 


(31) 


remains constant, as is easily verified using Eqs. ( |25p7t . 
Also, K is subject to the inequality K > ^, sometimes 
known as the Schrodinger-Robinson inequality [Refill 
Eq. (10.57)]. It is stronger than (and includes) Heisen¬ 
berg’s uncertainty relation A^Ap > H/2. An alternative 


proof of the the inequality will emerge in Section VHI 
From Eqs. (28) and (29) follows the relation 


= A^+K^ 


(32) 


and therefore e > K. Also e = AT if and only if A = 0 
which implie s no oscillations in the second moments. 
From Eq. (31), A^Ap = + A^ and A^p oscil¬ 


lates between —HA and HA. Therefore the ‘uncertainty 
product’ AxAp oscillates between HK and he. Thus 


e > A„Ap/fi > A > i. (33) 

A physical manifestation of K is that it is the minimum 
value that A^Ap/h takes as it oscillates. 
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Returning to the question of how to calculate the con¬ 
stants e, A and to in terms of the initial moments, e and 
K can be found from Eqs. (301 and (311. Then A follows 
from Eq. (32) while Eqs. (281 and (291 give (taking t = 0) 


Asm2ujto =—Axp/h (34) 

Acos2a;to = ^(^Ap-^A^). (35) 


VIII. CHANGING THE SCALE OE A STATE 


We saw that the evolution of any state can be reduced 
to that of a state that remains at rest at the origin. This 
transformation involved only a spatial translation and a 
change in phase that depends only on the first moments 
of the state. Now we show that the evolution of these 
states can be further reduced to states whose second mo¬ 
ments are constant, which we refer to as stable states. 


Stable states. Erom Eqs. ( 28p^ , having constant sec¬ 
ond moments requires A = 0. This implies that: K = e, 
A^p = 0, Al = A2 = h^K/a\ A,Ap = HK 

and these values will be maintained as the state evolves. 
Any state with Axp = 0 (which includes any with a real 
wave function) is stable if and only if A^fa = aApjh. 
Although the second moments of stable states do not 
change with time, the wave functions do change to some 
extent. Energy eigenfunctions are stable, but their phase 
does change with time. If if = | the stable state has 
e = \ and therefore must be the ground state Gaussian. 


Changing scale. In the case of a free particle, changing 
scale is a simple matter: if 'ip{xA) satisfies Schrodinger’s 
equation = —(^^/2to) 9^^ then so does s^t). 

Thus if we change the spatial scale, we still have a so¬ 
lution (with a different time scale). Eor the oscillator, 
we can find a solution with a different spatial scale, but 
this generally requires an added phase (that varies with 
X and t) and a change the time scale, now to a nonlinear 
function of the original time. 


Any state can be transformed to a stable form. 

Eor a state "(pix) with A^p = 0, which includes any with 
a real wave function (as shown in Appendix D), we can 
find the scale s such that ip^sx) has lower energy than 
any other scale, and this lowest-energy state is stable. An 
example is worked out for the triangular wave function 
in Example 2 of Section IX Otherwise, when ipix) has 


Axp ^ 0 we can find a phase change that reduces A 


xp 


to zero and then change the scale. The details are in 
Appendix C and the result is that 


(j){x) = exp(—^1x^/6^) V'(sa^), 


(36) 


is stable if s = Ax/a'/K and = a^hK/Axp, where 
Ax, Axp and K refer to ip{x). In fact K is unchanged by 
this transformation. 

Every stable state has e = K and its energy cannot be 
less than ^hio. Therefore every state must have A > I, 
which proves the strong uncertainty inequality. 


Evolution of a state from its stable form. The in¬ 
verse of Eq. p^, expressing ipix) in terms of 4>{x), is 
^p{x) = exp(|i^/6^s^)(/)(a;/s), and this is the starting 
point to find the evolution ip{x,t) of '4>{x) from the evo¬ 
lution (p{x,t) of (p{x). The time t in 4>(x,t) must be dis¬ 
torted by replacing t by some function r(t). Thus, if we 
insert 


“Vl) 


(37) 


_ 1 /o 

into Schrodinger’s equation (where the factor A^, ' is 
clearly required to preserve normalization), we find that 
ip{x,t) satisfies Schrodinger’s equation if (j){x,t) does, 
provided that dr/dt = Ka’^jA'l and A^,, Axp evolve as 
in Eqns.( 25p7| . [Similarly we can find the evolution of 
exp{—ixy2a^'ip{x) if we know that of any ipix).] 

The distorted time T(f) is found by integrating 
dr/dt = Ka^/A/x, inserting A^ from Eq. (29), to give 


1 e -\- A 

r(t) = — arctan [——— tana;(t — to)]- (38) 

u) K 

When A = 0 then e = K and hence r(t) = t — tg. An 
example of r(t) is shown in Eig. 3. As required by the 
periodicity properties, t equals t —to at multiples of T/4. 
[As an alternative, Eq. (38) is equivalent to exp2*a;r = 
(1 -I- tz)/{l — iz), where z =[[e + A)/K] tana;(t — to).] 



FIG. 3. 

The distorted time r(t) compared with the actual time t 
elapsed since a time when = 0. In the case shown, 
K — 1, A = 1 and therefore e = \/2. 


The evolution of a scaled stable state (/{sx, 0) from 
4>[x, t) is given by Eq. (371 initially taking Axp = 0. Then 


Axp = tiAsin2a;t and^A^ = a^(e — Acos2cut), so that 
ip{x, 0) = (j){sx, 0), where = K/{e — A) = {e + A)/K. 

Squeezed states result from changing the scale of the 
ground state. They have K — ^ and can be defined by 
the amplitude A of the moments. The evolution of the 
ground state is (j){x,t) = exp(—— ^tcot) and so, 
from Eq. (37), the evolution of the squeezed state is 


V’(*A) = ^M{rAxp - - l^xT{t)], (39) 
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where = a{e — Acos2wi)^/^ with Aq = (Aa;)t=Oj 
Aa;p = hAsin2ujt, u}T{t) = arctan[2(e + A) tanwt] and 
e = {A^ + This is easily extended to the evolution 

of excited energy eigenstates after a change of scale. 

For a study of the evolution of squeezed states using 
the propagator, see Ref. El Squeezed states are impor¬ 
tant in quantum optics, where the mathematics of the 
quantized oscillations of the electromagnetic field mirror 
that of the harmonic oscillator.l^ 


R:2 = 3/io, e=i(A2/10-p3/A2) and i(AVlO-3/A2), 
where A = a/a is the only non-trivial parameter for this 
wave function. Then the minimum energy, with e = K = 
(3/10)^/^ « 0.5477 occurs when A = 0 and requires that 
A = 30^/4 « 2.34. This is the stable triangular wave 
function. It is the wave function most similar (in this 
sense) to the ground state of the oscillator, and its second 
moments will remain constant. These two wave functions 
are compared in Fig. 4. 


Why change scale? Before carrying out a numerical 
calculation of the evolution of a wave function ip^x) (such 
as the triangular wave function in Example 2) it is useful 
to first transform the state to a stable form. Then the 
results can easily be applied to all the wave functions of 
the form 'ip^sx); all that is required is a relabelling of x 
and t and an easily calculated phase factor. 

More generally, calculations of physical processes often 
use a basis of states. The scale of the system may change 
significantly with time, and then it may give considerable 
advantage to allow the basis to change in scale to match 
the changing scale of the system. 


IX. EXAMPLES OF EVOLUTION 



FIG. 4. 

The stable (i.e. minimum-energy) triangular wave function 
(a « 2.34a) compared with the ground-state wave function. 
Both wave functions are normalized to dx = 1. 


Example 1 is a superposition of two ground-state Gaus¬ 
sian packets displaced from the origin, as shown in Figs. 1 
and 2. We will discuss how to calculate the evolution 
of the wave function and its first and second moments. 
The Gaussians start at rest, one centered on a; = oi 
while the other, with amplitude b relative to the first, 
centered on a; = 02 . Thus the initial wave function is 
tpix) = ipi{x) + btp2ix), where — a^) and 

'>p2{x) = <)>o(a; — 02 ), where (j)o{x) = exp{—x‘^/ 2 a^). The 
trajectory of the centroids will be Xk(t) = at cos{ujt) and 
Pk{t) = ttk cos{ujt) for fc = 1, 2. From Eq. (151, ipkix) will 
evolve to 


tpk{x,t) = exp[i9k{x,t) - 7^{x - Xk{t) f], (40) 

Za 


where 0k{x,t) = Pk{t){x — ^Xk{t))/h — Then 

ip(x,t) = ' 4 ’ii^A)+bip 2 ix,t) is the evolution of the super¬ 
position ijj{x) and is depicted in Eig. 1. The integrals for 
the initial second moments can be carried out exactly and 
the results contain only powers and exponentials (but are 
best done with a computer algebra/calculus program). 
Once the initial moments are calculated, Eqs. (281 and 
(291 give the evolution of the moments as shown in Eig. 2. 


Example 2. Another example is the triangular wave 
function, with ipi^) = for Is^l ^ o- In contrast to 

Example 1, the moments are easily calculated and reveal 
a lot about the evolution, but the exact evolution is not 
easily available. The wave function is real and it follows 
(see Appendix D) that initially {x) = {p) = {{xp+px))=0. 
It is easy to calculate that (x^) =a^/10 and {p'^)=3h^/a'^. 
Thus we have the initial values A^ = a^/10, A^ = ‘iP? jd^ 
and Aa;p = 0, from which we can calculate the constants 


The Fourier transform of the triangular wave function 
can also be calculated exactly and gives 


V^(x,T/4) 


ax"^ 


(41) 


For the minimum-energy triangular wave function, the 
absolute value of this Fourier transform does not differ 
greatly from the original triangle (as can be seen from the 
top graph in Fig. 5) and since Aj, is constant we might 
expect that the whole evolution remains close to the tri¬ 
angle. To be more specific, we can express the normal¬ 
ized initial state ')/' as a linear combination of the ground 
state ipo and a part y orthogonal to ipo- '0 = cq^o + 
with (V'l^) = 1, (0o|0o) = 1, (xlx) = 1 and (ipolx) = 0. 
Finding cq and requires only a numerical integration 
of '000 and the result is that Cg « 0.9953 and there¬ 
fore « 0.0047. Thus the non-Gaussian part of 0 con¬ 
tributes less than 1/200 to the probability distribution 
and this will be maintained throughout the evolution. 
The Gaussian part remains unchanged (apart from its 
phase) and all the deviations from the Gaussian come 
from this small x- 

This is borne out by a calculation of the evolution, first 
expanding the triangle as a sum of oscillator eigenfunc¬ 
tions (by numerical integration) and then including the 
time-dependence of the eigenfunctions. Figure 5 shows 
the result for |0(x,<)| over the first quarter period. 

For a wider initial triangle, that is a > the 

result is similar except that A^, decreases over the first 
quarter period. The evolving wave function stays close 
to an evolving Gaussian, which is now a squeezed state 
rather than a coherent state. 
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FIG. 5. 

The evolution of the minimum-energy triangular wave 
function over a quarter period. The time interval between 
each plot is r/32 and \ip\ is shown (solid curve) while the 
original triangular wave function is shown (dashed) at each 
time for comparison. 


X. CONCLUSIONS 


be evolved by first applying a simple transformation that 
depends only on (x) and (p). This greatly simplifies the 
practical problem of evolution of an oscillator. 

Without calculating the full evolution, a great deal 
about the behaviour of wave functions can be obtained 
from the three second moments A^, and A^p, which 
form a closed set: given their initial values we can find 
their evolution, independent of any other properties of 
the wave function. The energy can be separated into 
the energy Ec of the centroid (the energy of the equiva¬ 
lent classical oscillator) which depends only on the first 
moments, and the quantal energy Eq, which is the re¬ 
maining part of (H) = Ec + Eg. The quantal energy 
can be expressed in terms of second moments, Eq = 
Apl2m + and Eg > ^huj. Both Ec and Eg re¬ 

main constant. The evolution is also constrained by the 
constancy of the combination K'^ = (A^A^ — Alp)jh'^. 
This constant K is subject to the inequality K > ^, 
which contains (but is stronger than) the traditional 
Heisenberg uncertainty relation. The extreme case of 
K = ^ implies that the wave function is Gaussian, and 
is known as a squeezed state. Apart from the motion 
of the centroid and the choice of the origin of the time, 
squeezed states are characterized by a single parameter 
that defines the extent of the oscillations of the second 
moments. The case with no such oscillations is the coher¬ 
ent state, equivalent to the ground state of the oscillator. 


In principle, the evolution of oscillator wave functions 
is provided by the known propagator, but only a small set 
of initial wave functions allow the integral over the prop¬ 
agator to give a simple result and even then the main 
attributes of the evolution may be hidden in the math¬ 
ematical formulae. It is useful, therefore, to investigate 
some general features of the evolution. 

One simple feature is the periodicity. After one quar¬ 
ter of the classical period T = 27r/a;, any wave func¬ 
tion evolves to one that can be simply expressed in 
terms of the Fourier transform of the original wave 
function. A simpler result applies over half a period: 

^T) = —iiIj{—x,0). Applying this result twice gives 
tpix, T) = —4’{x, 0). Thus the true period of the quantum 
oscillator is twice that of the classical oscillator; evolution 
over 2T returns the wave function exactly to its original. 
These results imply that, for any initial wave function, 
the evolution over any time interval can be found from 
that through the first half period. If the initial wave 
function is real, only the first quarter period is needed. 

The simplest fact about the continuous evolution is 
Ehrenfest’s result that the centroid strictly follows clas¬ 
sical evolution; that is, sinusoidal oscillation that can be 
easily calculated from the initial values. The evolution of 
the whole wave function relative to the centroid, is essen¬ 
tially independent of that of the centroid, and therefore 
we need consider only states that remain centered on 
the origin (the center of force of the oscillator); that is, 
states with (x) = (p) =0. Any other wave function can 


For some states, here called stable states, the second 
moments do not change with time. If we have the evo¬ 
lution of a stable state tjj{x), then the evolution of any 
state '0(sa;) can easily be written down. Also, from any 
state ip{x) one can find a stable state through a transfor¬ 
mation that involves only its second moments, and the 
evolution of t/^ix) can be easily found from that of its 
stable state. Thus the evolution of every oscillator wave 
function (for which the second moments exist) is simply 
related to that of a stable state at rest at the origin. 


The main focus here has been the basic principles; but 
the potential applications are spread widely in physics, 
because harmonic oscillations are ubiquitous. In practice, 
there will be other physical processes acting (such as ra¬ 
diation and environmental effects). These can be dealt 
with using a basis of energy eigenfunctions of the oscil¬ 
lator, but in some cases it is more effective to shift and 
scale the basis to match the changes in the object being 
studied. The evolution of such displaced and sque ezed 
eigenstates was discussed in Sections 0 and dvml). 


In two or three dimensions all the results above can 
be applied in each dimension since the oscillator Hamil¬ 
tonian can be separated as a sum of commuting terms. 
Thus, for example, in the case of an isotropic oscillator in 
three dimensions, the centroid of every state will travel 
along an elliptic orbit in a plane and the evolution of 
the wave function can be reduced to that of a state that 
remains centered on the origin. 
































Appendix A: The propagator from an invariant 
operator 

The propagator K{x,x',t), as defined by Eq. Q, is the 
solution of Schrddinger’s equation such that K (a;, x',t) —> 
6{x — x') as t —>■ 0 and has the symmetry property 
K{x', x, —t) = K*{x, x’, t)P 
If we can find an operator u such that u —^ a; as t —>■ 0, 
then the solution of u'lp = x'lp^ will have "0 oc d{x — x') 
as t —?► 0, because S(x — x') is the eigenfunction of x with 
eigenvalue x'. If ii is also an invariant OMrator, that is 
Dtu = 0 with Dt defined in Eq. (101, theiP "0 will satisfy 


Schrodinger’s equation HiIj = ihdttp apart from a possible 
time-dependent factor, and will give the a;-dependence of 
the propagator. 

In the case of the oscillator, we can find such an oper¬ 
ator u as a linear combination of x and p. Thus if u = 
P{t)x — X{t)p, then DtU = {P + muj'^X)x— {X — P/m)p, 
so u will be invariant if X{t),P(t) satisfy the classical 
equations of motion. Therefore, we take 


u = cos uitx — (muj) ^ sin tot p. 


(AI) 


Then the solution of uijj = x'ip is 


exp[ 


21isina;t 


{cosiutx^ — 2x'x)], (A2) 


and this gives the cc-dependence of the propagator. The 
symmetry between x and x' then requires 


K{x,x\t) = / 2 (t)exp[ 


imu! 
2H sin Lot 


(cos Ljt{x'^ +x''^) — 2x'x )], 


and inserting this into Schrddinger’s equation shows that 
K{x,x',t) will be a solution provided that dtf 2 {t) = 
—w cotwt/ 2 (t) and therefore / 2 (t) = c/psinut. All 
that remains to be determined is the constant c and 
this can be found by requiring that, for small t, 
X(x, x', t) exp (—dx should approach unity as 

b —>■ oo. This gives / 2 (t) = y/rnuo/2'Kih sin cot. Making 
the phase precise, the propagator is therefore 


K{x, x', t) = 


e-x7r/4^-l 

1^2711 sinwtl 


exp[i 


{x^ + x'^) cos Lot — 2xx' 
20^^ sin Lot 


(A4) 

agreeing with Eq. Q . As discussed in Section |III[ the 
phase given here is valid only for 0 < Lot < ir because of 
the singularity when sin cat = 0. 

Resolving the ambiguity in phase. This analysis 
shows that a representation of the (5-function is 


5{x) = lim 


g— 27r/4 


exp 


IX 


e^o+ aV2jLx 2eQ(2 


(A5) 


so that, fo r sm all oot, K{x,x',t) —>■ S{x — x'). For oot = 
TT — e, Eq. (A4) becomes 


K{x,x',t) « ——— exp 


e —i{x + x') 


,/^2 


(A6) 


aV^e " 2ea2 
e“*’’’/^(5(a; -I- x') as e —>■ 0, (A7) 


where we have used the complex conjugate of Eq. (A5) 
for the (5-function. Hence, ip{x,t) « a;, 0), as 

required. Now if we pass through the singularity to time 
Lot = TT + e, Eq. (A4) gives K{x,x',t) = 6{x + x'), which 
is incorrect as it involves a sudden jump in phase. To 
get the correct result, we must multiply K{x,x',t) by 
g-iir /2 _ This extra phase is required each time cot 


passes through a multiple of tt. Therefore Eq. (A4) re¬ 
quires an extra factor of (—*)^, where k is the integer part 
of Lotjir. This phase change can be better understood 
through Maslov theory and phase space geometry.® 


Appendix B: The Fourier transform of 
eigenfunctions and evolution over r/4 

The standard lowering operator for the oscillator is a = 
{x/a -\- iap/h)/y/2 and the ground state ipo has S'l/'o = 0 
with the wave function 'ipo{x) = exp{—x^/2a'^). From 
this ground state, we can generate the infinite sequence 
of excited eigenstates by repeatedly applying the raising 
operator d^. Thus ipn+i = c„d^(/5„, where c„ is some real 
number that could be taken to preserve normalisation, 
but is not relevant here. 

We will now find a relation between 'ipn{x) and its 
Fourier transform (or the momentum wave function). 
This relation arises from the symmetry between posi¬ 
tion and momentum in the oscillator, and it helps to 
use dimensionless variables ^ = x/a and p = ap/h 
to make this symmetry more explicit. [For example, 
(A3) Schrodinger’s equation becomes -I- ?)'*/’ = 


Write ip{x) = '^{x/a) and define its Fourier transform as 


$(p) =p-^ = 


1 


y/2TT J- 


/ OO 

(Bl) 

-OO 


Then the ground state = exp(—has the 

Fourier transform $o(p) = exp(—fp^). Now suppose 
we have some function 'l'(^) and examine the Fourier 
transform of d^'k(^). We have = (^ — d^)/\/2 and, 
integrating by parts, 

/ OO nOO 

-oo J —OO 

/ OO 

(B2) 

-OO 

That is, 

j:-(d1’4') =-7d^(J-4'). (B3) 


We showed that = 'I'o(p); then Eq. (B3) gives 

•^ 4 ^ 1(0 = —*4'i(p). [If we take the first excited eigen¬ 
function to be 'ki(^) = d^4>o(0/-\/2 = ^exp(—4^^), then 
its Fourier transform is $i(p) = —jpexp(—4p^).] This 
process can be continued, and an extra —i is factored in 
at each step. Thus 

■^4'„(0 = (-*)”4/„(p). 


(B4) 
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More formally, if Eq. ( |B4[ ) holds for some n, then Eq. (B3) 
gives = (—and therefore 


We now apply this result to the evolution of any wave 
function over a quarter period. As for Eq. ([^, expand the 
wave function as a sum over the energy eigenfunctions: 
4'(^) = J2^=o C'n^n(0- Then the evolved wave function 
at time t is 


E C50 ^ 

exp[-iwt{n + -)]C'„4'„(C). (B5) 

n—0 Z 


If t = T/4 = 7r/2a;, the factor exp[—ia;t(n + i)] becomes 
exp(—z7r/4)(—i)" and therefore 


It follows that K is unchanged and that we can make 
^xp = 0 by taking = hA'^/s'^Axp. Then e becomes 
lis^a^XyAl + Al/s^ a^) and takes its minimum value 
e = K with s = A^ja'/K. Hence = Ha^K/A^p and 


This wave function (^(x) is the stable equivalent of il’ix)- 


T/4) = (B6) 


Appendix D: Properties of real wave functions 


But the Fourier transform of 4'(^) is 4?(p) = 

= E“oC'n-^'I^n(0 = 

Thus 'I'(^,T/4) is, apart from the phase shift of 7r/4, the 
same function of ^ as is a function of p, agreeing 

with the result from the propagator in Section m 

Appendix C: Transforming a state to stable form 

We show that changing the scale of a state cannot yield 
a stable state unless A^p = 0; but after a suitable phase 
change, a stable form can always be reached. 

If the moments of ijj{x) are A^,, A^p, Ap, and we 
change the scale of tpix) to give x(a;) = ip{sx), where 
s is real, then it follows from the integrals for the expec¬ 
tation values that the moments of x{x) are = s~^Ax 
and Ap = s^Ap, while Aa,p, A^Ap and K do not change. 
Then e becomes is“^A^/a^ -|- ^s'^a^Ap/h? and this has 
a minimum value of AxAp/h, = + A^pjfi^y-/'^ at 

= hA^/a^Ap. The lowest energy depends on A^-p and 
only if Aa;p = 0 can we reach e = K, which implies A = 0. 

We can change A^p by a change in phase. Direct calcu¬ 
lation gives the moments of (/)(x) = exp{—ix'^/2b'^)il;{sx): 

Al = Al/s\ (Cl) 

A,p = A,p-hAl/sH\ (C2) 

Al = s^Al - 2H A,p/&2 + h^Alls%^. (C3) 


Many textbook examples have initial wave functions 
that are real, and therefore have special properties. 
These properties are ephemeral because every real wave 
function will immediately become complex (due to the 
term ihdtip in Schrodinger’s equation). However some 
special properties of real wave functions are useful in the 
present context. 

First we show that (p) = 0 when ip is real. We use 
the Hermitean property of p = —ihdx- J(p*pxdx = 
J{p(p)*xdx for any two wave functions (p and y. Hence 
Jip*pipdx = J(j>ip)*ipdx = J {idxip)ip dx = — JIp pip dx, 
since ip* = ip. Similarly, A^-p = 0: Jippxipdx = 
f {pip)*xip dx = — f {pip){xip) dx = — Jip xpip dx. 

Unless {x) = 0 also (in which case the centroid sits 
permanently at the origin), reality of the wave function 
implies that the centroid is instantaneously at rest at 
the extreme of its travel while A^, and Ap also take their 
extreme values. Thus the phases of the oscillations of the 
centroid are aligned with those of the second moments. 

We can add momentum to any state (p using wave 
functions of the form ip = exp(ip'x/h)(p. In this case, 
{x)^ = {x)rp while pip = p'lp + e^'P ^^^p(p and there¬ 
fore (p)^ = p' + {p)cj,. Also (xp)^ = p'{x)p + {xp)<j, 
while (px) = (xp) — ih for either ip or (p, and therefore 
{px + xp).^ = {px + xp)^ -\- 2p'{x)^. If (p is real, it follows 
that {Axp')i^ — 0. 
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